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5 APPARATUS AND METHOD FOR IMPLEMENTING EFFICIENT ARITHMETIC 
CIRCUITS IN PROGRAMMABLE LOGIC DEVICES 



10 

This application claims priority to the provisional patent application entitled, 
"Apparatus and Method for Implementing Complex Arithmetic Circuits in Programmable 
Logic Devices," serial number 60/236,244, filed September 28, 2000. 

1 5 BRIEF DESCRIPTION OF THE INVENTION 

The present invention relates generally to programmable logic devices (PLDs). 
More specifically, the present invention relates to techniques for efficiently implementing 
complex circuits in PLDs. 

20 BACKGROUND OF THE INVENTION 

Many applications in digital signal processing require the use of complicated circuits 
for performing complex arithmetic calculations in real time. For example, in order to cancel 
echoes that occur as a result of hybrid coupling within a telephone network, adaptive Finite 
Impulse Response (FIR) filters have been used. Such FIR filters typically include 
25 complicated arithmetic circuits for performing reciprocal calculations, and complicated ^i- 
Law/A-Law expander circuits for expanding |Li-Law/A-Law encoded signals. 

PLDs are the devices of choice in implementing these complicated digital signal 
processors. For many PLDs, the basic building blocks are the Logic Elements (LEs) or 

30 Logic Array Blocks (LABs) that include logic circuits that may be programmed to perform 
specific logic operations. Due to their complexity, digital signal processors typically require 
the use of many LEs. Each PLD has a limited number of LEs. Therefore, after these 
arithmetic circuits are implemented, few LEs are left for performing other fimctions. As a 
result, it is often difficult to implement an entire digital signal processing system with 

35 compHcated arithmetic circuits on a single PLD. 
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Accordingly, it would be desirable to provide improved techniques for efficiently 
implementing complex arithmetic circuits in programmable logic devices. 

SUMMARY OF THE DISCLOSURE 

5 Techniques for efficient implementation of complex arithmetic circuits in 

programmable logic devices are disclosed. In one aspect of the present invention, Look-Up 
Tables (LUTs) of a programmable logic device are used to store pre-calculated intermediate 
or final calculation values. A table look-up operation is performed in place of complex 
arithmetic operations, hi this way, at the expense of a few LUTs, many logic elements can 

10 be saved. This approach is particularly apphcable to circuits for calculating reciprocal 
values and circuits for performing a normalized Least-Mean-Squared (LMS) algorithm. 

According to one aspect of the invention, a circuit for calculating reciprocal values 
is implemented with LUT(s) having stored therein pre-calculated reciprocal values. The 

15 pre-calculated reciprocal values are determined based on the address locations at which they 
are stored, hi operation, an input value to the circuit is partitioned into a number of 
segments. Then, one of the segment is selected for indexing the LUT(s) and for retrieving 
an associated pre-calculated reciprocal value therefrom. The retrieved pre-calculated 
reciprocal value is then shifted, according to a position of the selected segment within the 

20 input value, to obtain an output value. In one embodiment of the invention, the selected 
segment is the first segment, starting from the decimal pouit, that contains at least one 
digital "high" value (e.g., a logic "1" ). Further, according to the present invention, the 
input value is a fraction, and the output of the circuit, which is an approximate reciprocal of 
the input value, is an integer. 

25 

In another aspect of the present invention, a circuit implementing a normahzed LMS 
algorithm has LUTs storing pre-calculated |i-Law or A-law expansion values. An input to 
the circuit, which is either pi-Law or A-Law encoded, is provided to the LUT as an index for 
retrieving an associated expansion value. The retrieved expansion value is then provided to 
30 circuits for performing other steps of the normahzed LMS algorithm. Accordingly, by using 
LUTs to store pre-calculated [i-Law or A-Law expansion values, such complicated 
arithmetic circuits can be implemented with a significantly smaller number of LEs. 



35 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of the invention, reference should be made to the 
following detailed description taken in conjunction with the accompanying drawings, in 
5 which: 

Figure 1 is a block diagram illustrating a programmable logic device on which 
embodiments of the present invention maybe practiced. 

Figure 2A is a block diagram illustrating a circuit for calculating reciprocal values in 
1 0 accordance with one embodiment of the present invention. 

Figure 2B is a block diagram illustrating a circuit for calculating reciprocal values in 
accordance with another embodiment of the present invention. 

J 5 Figure 3 is a flow chart diagram illustrating a method of calculating reciprocal 

values using an LUT-based reciprocal calculation circuit in accordance with an embodiment 
of the present invention. 

Figure 4A is a block diagram illustrating a Look-Up Table (LUT) for implementing 
20 a |i-Law expander in accordance with one embodiment of the present invention. 

Figure 4B is a block diagram illustrating a Look-Up Table (LUT) for implementing 
an A-law expander in accordance with one embodiment of the present invention. 

25 Figure 5 is a block diagram illustrating a Look-Up Table (LUT) in accordance with 

an embodiment of the present invention. 

Like reference numerals refer to corresponding parts throughout the drawings. 

30 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In the following detailed description, for purposes of explanation, numerous specific 
details are set forth in order to provide a thorough vmderstanding of the present invention. 
However, it will be apparent to one skilled in the art that the present invention may be 
practiced without these specific details. In other instances, well-known structures, devices, 
35 methodologies, etc., are not described in detail in order to avoid obscuring aspects of the 
present invention. 
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A. General Description of a Programmable Logic Device 

Figure 1 is a block diagram illustrating a programmable logic device 100 (e.g., 
Altera FLEX™ lOK device). The device 100 includes logic arrays blocks 120, embedded 
array blocks 140, an interconnect grid 160, and input/output elements 180. Interconnect 

5 grid 160 includes sets of horizontal conductors 162 and sets of vertical conductors 164 that 
are configiired in a grid pattern. Logic array blocks 120 and embedded array blocks 140 are 
programmably coimectable to horizontal conductors 162 and vertical conductors 164 of 
interconnect grid 160. Interconnect grid 160 is also programmably connectable to 
input/output elements 180. By programming a connection between input/output elements 

10 180, logic array blocks 120 and embedded array blocks 140 through horizontal conductors 
162 and vertical conductors 164, the user may create desired logical functions. Particularly, 
PLD 100 maybe programmed to implement complex arithmetic circuits of the present 
invention. 

15 B. Implementation of Arithmetic Circuits using Look-Up Tables 

According to the present invention, complicated arithmetic circuits may be 
efficiently implemented within a programmable logic device (e.g., device 100) by using 
Look-Up Tables (LUTs) to store pre-calculated values. A table look-up operation can then 
be performed in place of complex arithmetic operations. In this way, at the expense of a 
20 few LUTs, many logic elements can be saved. 

As an example, it is desired to implement a circuit for calculating a complicated 
arithmetic function fii(x) in response to an input value x. Generally, comphcated logic 
circuits requiring a large number of logic elements are used. However, according to an 

25 embodiment of the present invention, values of the function fii(x) over a range of input 
values (x) are first determined. The pre-determined or pre-calculated values are then stored 
within a LUT of a programmable logic device. In one embodiment, the pre-determined 
values are stored at addresses corresponding to the input values (x). For example, the pre- 
determined value of fh(l) is stored at address 1, and the predetermined value of fh(255) is 

30 stored at address 255. The LUT can then be used for calculating the value of fh(x) for an 
input value (x). For example, when an input value k is presented to the look-up table, the 
value stored at address k, which is fli(k), will then be provided as an output. A block 
diagram of a LUT 500 and its contents according to this embodiment of the present 
invention is illustrated in Figure 5. The circuit for matching the input value to the 

35 appropriate address of the LUT 500 is well known in the art, and is not illustrated for 
purposes of clarity. 
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The present invention is particularly applicable to the implementation of 
complicated arithmetic circuits for adaptive digital voice filtering in programmable logic 
devices. In the following, exemplary implementations of circuits for carrying out 
normalized LMS (Least-Mean-Squared) algorithms are presented. It should be appreciated 
5 that, however, the present invention is not limited to such circuit implementations. 

C. Implementation of a Reciprocal Circuit in a Programmable Logic Device 

In various digital signal processing appUcations, efficient calculations of reciprocal 
values are often required. A brief description of a fast LUT-based algorithm for computing 
IQ a reciprocal value given a fraction with arbitrary bit width follows. 

Suppose P„ = .(an., an.2a^.3 .... aj ao) is an n-bit binary unsigned fraction. The decimal 
value of Pn can be expressed as follows: 

15 Pn = a„., *2-' + a„.2*2-^ + a„.3*2-^ .... + ao*2-" 

(a„.i * 2-^ + a„.2 * 2'^ + a„.3 * 2'^ + ... + a^.L * 2'^) + 

2-" (a „.(L.i) * 2-' + a „.(lh-2) * 2"^ + a „.(l.3) * 2"^ + ... + a * 2"^) + 

2 n-(2L+l) * 2 ' + a n-(2L+2) * 2 + E n-(L+3) * 2 ^ + .... + a „.3l * 2^) + .... 

2 --"^ (a„.(.L.i) * 2-' + a„.^^,2) * 2"^ + a„,^L.3) * 2"^ + ....+ ao * 2"^). 

20 

A certain symmetry can be seen from the above expression. This symmetry of the 
expression can be exploited. Particularly, all possible reciprocal values of (a * 2'' + a „.2 * 
2"^ + a „.3 * 2'^ + ... + a „.l * 2'^) can be pre-calculated and stored within a look-up table, 
using the co-efficients (a n.(L+i) a n.(L+2) ^ n-(L+3) — a ^-21) as addresses. The reciprocal value of 

25 2'^ (a „.(L+i) * 2-' + a „.(l+2) * 2"^ + a „.(l+3) * 2"^ + ... + a „.2l * 2"^) can be determined by 

retrieving the pre-calculated reciprocal value from the look-up table using the coefficients (a 
ii-(L+i) a n.(L+2) a n-(L+3) — a „.2l) as the address, and by multiplying the retrieved value by 2^. 
Similarly, the reciprocal value of 2 (a„.(„L+,) * + a„.(^L+2) * 2"' + a„.(n,L+3) * 2"^ + ....+ a^ * 
2'^) can be determined by retrieving the pre-calculated reciprocal value from the look-up 

3Q table using the coefficients (a^^^^L+i) an.(mL+2) an.(inL+3) ••• a^) as the address, and by multiplying 
the retrieved value by 2 

As an example, suppose the fraction P„ is divided into m equally wide fields each 
having 4 bits. To calculate the reciprocal of P„, m look-up tables each having 4-bit wide 
35 addresses (or indices) and each capable of storing sixteen reciprocal values can be used. 
The contents of one such look-up table are illusfrated below in Table 1 . 
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Table 1 



15 



20 



Address 


Content 


0000 


0 


0001 


16 


0010 


8 


0011 


16/3 


0100 


4 


0101 


16/5 


0110 


8/3 


0111 


16/7 


1000 


2 


1001 


16/9 


1010 


8/5 


1011 


16/11 


1100 


4/3 


1101 


16/13 


1110 


8/7 


nil 


16/15 



In operation, for a fraction .1011 (bin), its reciprocal value can be determined by 
looking up Table 1 to retrieve the value stored at address 101 1, which is 16/1 1, or 
approximately 1.23. For a fraction .00001001 (bin), its reciprocal value can be calculated 
by looking up Table 1 to retrieve the value stored at address 1001, which is 16/9, and 
multiplying 2^ to the retrieved value. In the present example, the reciprocal value of 
.00001001 (bin) can be easily determined to be 2"* * 16/9, or approximately 28.44. 

Reciprocal calculation can be sped up even further by storing the reciprocal values 
in the look-up table as binary integer values. Multiplication of binary integer values by 2" is 
equivalent to shifting the values to the left by n bits. Naturally, accuracy of the reciprocal 
calculation may be somewhat compromised. Nevertheless, in some appUcations where 
speed or simplicity is desired, such a trade-off may be preferred. 
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Note that the least significant segment may include fewer than L-bits (e.g., n-1 mod 
L bits). In that case, the look-up table with fewer entries will be used. It should also be 
noted that the present invention may also be applied to signed values. In that embodiment, 
the magnitude of a signed value is first determined. Then, the reciprocal of the magnitude is 
determined using techniques discussed above. Thereafter, if the original number is 
negative, the reciprocal value is two-complemented. 

Figure 2A is a block diagram illustrating a circuit 200 for calculating reciprocal 
values in accordance with one embodiment of the present invention. Circuit 200, according 
to the present invention, may be implemented using a programmable logic device (e.g., 
device 100). As illustrated, circuit 200 includes 2's complement circuits 210a-210b, a 
plurality of identical look up tables 220a-220n, a plurality of shifters 230a-230n, and an 
output selector circuit 240. 

In operation, when a number P,, is received, the 2's complement circuit 210a 
determines whether the number is negative and determines the magnitude | of the 
number. The magnitude jP^ | is then partitioned into m segments each being L-bit in width. 
The partitioned magnitude iP„ | 225, containing a plurahty of segments 235a-235n, is also 
illustrated in Figure 2A. The most significant bit of the partitioned magnitude |P„ | is the 
sign bit 234, which is not taken into account when the number is partitioned. Note that the 
least significant segment 235n may contain only (n-1 mod L) bits. 



The segments 235a-235n are then provided to look-up tables 220a-220n. According 
to the present embodiment, each LUT 220a-220n has stored therein a plurahty of pre- 
determined reciprocal values. In the embodiment illustrated in Figure 2A, the contents of 
the LUTs 220a-220n may be identical. Further, the reciprocal values are pre-calculated 
based on the addresses of the LUT at which the reciprocal values are stored. For example, 
address 1011 stores the reciprocal value of . 1 0 1 1 , and address 1110 stores the reciprocal 
value of . 1 1 1 0, etc. Thus, when a segment 235a of the magnitude IP^ | is presented to the 
LUT 220a, the LUT 220a retrieves the reciprocal value stored at the address corresponding 
to the segment 235a and provides the retrieved reicprocal value to the shifter 230a. 
Likewise, when a segment 235b is presented to the LUT 220b, the reciprocal value stored at 
the address corresponding to the segment 235b will be retrieved and provided to the shifter 
230b. It should be noted, however, that address 0000 does not store the reciprocal value of 
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.0000. Rather, address 0000 of the look-up table stores a value of zero, or a value that 
indicates a potential error. 

Each of the shifters 230a-230n shifts the retrieved reciprocal values according to the 
^ position of the segment vi^ithin !?„ |. For example, the shifter 230b is configured to shift the 
reciprocal value L bits to the left, and the shifter 230n is configured to shift the reciprocal 
value (n-1) * L bits to the left. 

The shifted values are then provided to the output selector circuit 240. In 
J Q accordance v^^ith the present embodiment, the output selector circuit 240 selects a first non- 
zero output from the shifters 230a-230n and ignores the rest of their outputs. In this 
embodiment, the shifter 230a, corresponding to the most significant bits of the input value, 
has the highest priority, and the shifter 230b, corresponding to the next most significant bits 
of the input value, has the second highest priority, etc. For example, if the output of shifter 
230a is non-zero, then the outputs of the shifters 230b-230n are ignored. As another 
example, if the output of shifter 230a is zero (indicating that the first segment 225a contains 
0000), but the output of the shifter 230b is non-zero, then the selective adder 240 selects the 
output of the shifter 230b. The selected output is then provided to 2's complement circuit 
210b to be converted into 2's complement form if the original niimber is negative. 

20 

In another embodiment of the present invention, instead of using multiple LUTs in 
parallel, a single LUT may be used to determine the reciprocal values. Figure 2B illustrates 
a circuit 205 for calculating reciprocal values in accordance with another embodiment of the 
present invention. Like circuit 200, circuit 205 is suitable for implementation by a PLD. 

25 

As shown, circuit 205 includes 2's complement circuits 210a-210b, an input segment 
selector circuit 250, a look-up table 260, and a shifter 270. The LUT 260 has stored therein 
a plurality of pre-determined reciprocal values that are pre-calculated based on the addresses 
of the LUT at which the reciprocal values are stored. For example, address 1011 stores the 
reciprocal value of . 1 0 1 1 , and address 1110 stores the reciprocal value of . 1 1 1 0, etc. 

In operation, when a nimiber P„ is received, the 2's complement circuit 210a 
determines whether the number is negative and determines the magnitude |Pj, | of the 
number. The magnitude [P„ j is then partitioned into m segments each being L-bit in width. 
The partitioned magnitude \P„ \ 225, containing a plurality of segments 235a-235n, is also 
illustrated in Figure 2A. The most significant bit of the partitioned magnitude |Pj, | is the 
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sign bit 234, which is not taken into account when the number is partitioned. Note that the 
least significant segment 23 5n may contain only (n-1 mod L) bits. 



The segments 235a-235n are then provided to the input segment selector circuit 250. 
The input segment selector circuit 250 selects the first non-zero segment of the magnitude 
|P„ I 225 (e.g., the non-zero segment that is closest to the sign bit 234) to be provided to the 
LUT 260. In response, the LUT 260 retrieves the reciprocal value stored at the address 
corresponding to the first non-zero segment and provides the retrieved value to the shifter 
circuit 270. The shifter 270 then shifts the retrieved reciprocal value according to a shift 
control signal provided by the input segment selector circuit 250. The output of the shifter 
270 is then provided to 2's complement circuit 210b to be converted into 2's complement 
form if the original number is negative. 



Figure 3 is a flow diagram illustrating a method of determining reciprocal values for 
an input data using a LUT-based reciprocal calculation circuit (e.g., circuit 205) in 
accordance with one embodiment of the present invention. As shown in Figure 3, at step 
310, reciprocal values are stored in a look-up table. According to an embodiment of the 
invention, the reciprocal values may be stored at addresses of the look-up table that 
correspond to the reciprocal values. For instance, the reciprocal value of .0010 may be 
stored at address 0010 of the LUT. It should be noted that, in that embodiment, an ERROR 
value or a zero value maybe stored at address 0000 of the LUT. 



At step 320, an input data is received. Then, at step 330, the input data is partitioned 
into a number of segments each having a bit-width corresponding to the bit-width of the 
addresses of the look-up table. For instance, if the look-up table has 8-bit wide addresses 
and the input data is 32-bit wide, the input data is partitioned into four 8-bit wide segments. 

At step 340, one of the segments is selected and provided to the look-up table to 
determine if there is a match and to retrieve a reciprocal value whose address matches the 
selected segment. According to the present embodiment, the first segment containing a 
non-zero value may be selected. Note that other selection criteria may also be used. 

Then, at step 350, the retrieved reciprocal value is shifted according to a position of 
the selected segment in relation to the input data to generate an approximate reciprocal 
value for the input data. For instance, if the second 8-bit segment of a 32-bit input fraction 
is selected, then the retrieved reciprocal value is shifted eight bits to the left. If the third 8- 
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bit segment of the 32-bit input fraction is selected, then the retrieved reciprocal value is 
shifted sixteen bits to the left. The result of the shifting will be an approximate value of the 
reciprocal of the 32-bit input fraction. 

D. Implementation of a |ii-Law/A-Law Expander Circuit 

The normalized LMS algorithm is a well-known algorithm in the field of 
telecommunications. Specifically, the normaUzed LMS algorithm can be used for canceling 
echoes that occur as a result of hybrid coupling within a telephone network. Generally, the 
normalized LMS algorithm can be represented by the following equations: 
P(j) - (1-a) * P a-1) + a * xG) * x(i) 
hiO+l) = hjG) + {[I ' eO)/PCj)) • x(j-lX where i = 0, N-1. 



For n-Law encoded data RinO"), the linear input data x(j) can be represented by the 
expression: |i-Law(Rin(j)), where Rin(j) is the |j,-Law 8-bit compressed value of the input 
data x(j). For A-Law encoded data Rin(j), the linear input data x(j) can be represented by 
the expression: A-Law(JUn(j)), where Rin(j) is the A-Law 8-bit compressed value of the 
input data x(j). 



Due to its complexity, the normalized LMS algorithm, if implemented as a circuit 
within a programmable logic device, would require more than six hundred logic elements. 
Logic elements, however, are critical resources in a programmable logic device. Thus, it is 
desirable to implement a circuit for performing LMS algorithms that do not require as many 
logic elements. The present invention presents a new approach that is significantly more 
efficient than prior art approaches. 

As shown above, the first step in the normalized LMS algorithm is to square the far 
end linear input data x(j) for calculation of the input signal power P(j). According to the 
present invention, instead of first calculating the linear value of x(j) = |j,-Law(Rin(j)) or A- 
Law(Rin(j)), and then calculating the value of x(j) * x (j), pre-calculated values for 
Law(Rin(j)) * |J.-Law(Rin(j))] or [A-Law(Rin(j)) * A-Law(Rin(j))] are stored in a single 
look-up table. In particular, the pre-calculated values are stored in the look-up table using 
the values (Rin(j)) as storage addresses. For example, pre-calculated value for [^i-Law(255) 
* !X-Law(255)] may be stored at address 255 of the look-up table. Then, during the LMS 
calculation, the input value Rin(j) is then used as an index to look-up the corresponding [^- 
Law(Rin(j)) * n-Law(Rin(j))] or [A-Law(Rin(j)) * A-Law(RinCj))] value. In this way, a 
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significant number of logic elements can be saved at the expense of only one embedded 
array block, which is used for implementing the look-up table. 

Figure 4 A is a block diagram illustrating a look-up table 410 for implementing a 
^ Law expander in accordance with one embodiment of the present invention. As shown, the 
look-up table 410 includes pre-calculated values of [|j.-Law(Rin(j)) * |Li-Law(Rin(j))] for all 
possible input values RinO'). It should be noted that Rin(j) is an 8-bit value. Thus, in the 
present embodiment, 256 pre-calculated jx-Law expansion values are stored within the look- 
up table 410. 

10 

Figure 4B is a block diagram illustrating another look-up table 420 for 
implementing an A-law expander in accordance with one embodiment of the present 
invention. The look-up table 420 includes pre-calculated values of [A-Law(Rin(j)) * A- 
Law(Rin(j))] for all possible input values Rin(j). Similar to table 410, the table 420 
^ ^ includes 256 pre-calculated A-Law expansion values. 

The present invention, the use of embedded memory for efficient implementation of 
complex arithmetic circuits, has thus been disclosed. The foregoing descriptions of specific 
embodiments of the present invention are presented for purposes of illustration and 
description. They are not intended to be exhaustive or to Hmit the invention to the precise 
forms disclosed, obviously many modifications and variations are possible in view of the 
above teachings. The embodiments were chosen and described in order to best explain the 
principles of the invention and its practical apphcations, to thereby enable others skilled in 
the art to best utilize the invention and various embodiments with various modifications as 
2^ are suited to the particular use contemplated. It is intended that the scope of the invention 
be defined by the following claims and their equivalents. 
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